Background
This notebook is for preparing input files for pyclone-vi. The same
files with a minor modification will be used as an input for FastClone.
Set up
suppressPackageStartupMessages({
library(tidyverse)
library(readxl)
library(cDriver) # Calculate CCF, https://github.com/hanasusak/cDriver/
})
Load and process data
# Read maf
maf_df <- readr::read_tsv(maf_file, guess_max = 100000, show_col_types = FALSE) %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, cg_id, tumor_descriptor,
sample_id, aliquot_id, Chromosome,
Start_Position, Reference_Allele, Tumor_Seq_Allele1, Tumor_Seq_Allele2,
t_ref_count, t_alt_count, tmb,
tumor_fraction, tumor_ploidy, VAF, Hugo_Symbol,
mutation_count) %>%
# Remove hypermutants
filter(!tmb >= 10,
# There are alterations with high number of read counts.
# We will exclude those with >1000 for now.
!t_alt_count >= 1000,
!t_ref_count >= 1000) %>%
# Add `normal_cn`: Total copy number of segment in healthy tissue.
# For autosome this will be two and male sex chromosomes one.
# See: https://github.com/Roth-Lab/pyclone-vi
mutate(normal_cn = case_when(grepl("chrY", Chromosome) ~ "1",
TRUE ~ "2"))
# Calculate CCF
maf_df <- CCF(maf_df)
# Read patient list
genomic_paired_df <- readr::read_tsv(genomic_paired_file, guess_max = 100000, show_col_types = FALSE) %>%
select(!c(cancer_group, experimental_strategy)) %>%
left_join(maf_df, by = c("Kids_First_Participant_ID")) %>% # create unique identifiers
dplyr::mutate(mutation_id = paste(Kids_First_Participant_ID, tumor_descriptor, Kids_First_Biospecimen_ID, Chromosome, Start_Position, Reference_Allele, Tumor_Seq_Allele2, sep = ":"),
cg_id_kids = paste(cg_id, Kids_First_Participant_ID, sep = "_"),
cg_id_kids = str_replace(cg_id_kids, "/|-", "_"),
cg_id_kids = str_replace_all(cg_id_kids, " ", "_"),
cg_id = str_replace(cg_id, "/|-", "_"),
cg_id = str_replace_all(cg_id, " ", "_"),
sample_id = paste(tumor_descriptor, Kids_First_Biospecimen_ID, sep = ":"))
# Let's count #specimens per sample
# We will remove any samples with less than 2 specimens as phylogenies require at least 3 taxa
kids_specimens_n_df <- genomic_paired_df %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, tumor_descriptor) %>%
unique() %>%
dplyr::count(Kids_First_Participant_ID) %>%
dplyr::mutate(kids_specimens_n = glue::glue("{Kids_First_Participant_ID} (N={n})")) %>%
dplyr::rename(kids_specimens_number = n) %>%
filter(!kids_specimens_number <= 2) %>%
left_join(genomic_paired_df, by = c("Kids_First_Participant_ID")) %>%
#left_join(maf_df) %>%
filter(!is.na(Chromosome))
# Let's confirm and add the number of timepoints per sample
# In case we lost any from filtering hypermutants and high reads/alteration
timepoints_n_df <- kids_specimens_n_df %>%
select(Kids_First_Participant_ID, tumor_descriptor) %>%
unique() %>%
dplyr::count(Kids_First_Participant_ID) %>%
dplyr::mutate(kids_timepoints_n = glue::glue("{Kids_First_Participant_ID} (N={n})")) %>%
dplyr::rename(kids_timepoints_number = n) %>%
filter(!kids_timepoints_number <= 1)
df <- timepoints_n_df %>%
left_join(kids_specimens_n_df, by = c("Kids_First_Participant_ID")) %>%
write_tsv(file.path(results_dir, "samples_eligible_for_phylogeny.tsv"))
# List with samples eligible for phylogeny
# I added the information about to use for phylogenetic inferences
list_df <- df %>%
select(Kids_First_Participant_ID) %>%
unique() %>%
mutate(somatic_germline_phylogeny = case_when(grepl("PT_KZ56XHJT|PT_KTRJ8TFY", Kids_First_Participant_ID) ~ "yes",
TRUE ~ "not_yet")) %>%
write_tsv(file.path(results_dir, "samples_eligible_for_phylogeny_list.tsv"))
Add Nautilus location for Deceased specimens
nautilus_dec_df <- read_excel(nautilus_dec_file) %>%
right_join(df, by = c("sample_id", "aliquot_id", "tumor_descriptor")) %>%
write_tsv(file.path(results_dir, "nautilus_dec.tsv"))
list_df <- nautilus_dec_df %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, `Note field from Nautilus of initial parent`) %>%
unique() %>%
write_tsv(file.path(results_dir, "nautilus_dec_list.tsv"))
Other, maybe to delete
# Make list with samples with >2 biospecimens at Deceased timepoint
dec_n_df <- df %>%
filter(tumor_descriptor == "Deceased") %>%
select(Kids_First_Participant_ID, tumor_descriptor, Kids_First_Biospecimen_ID) %>%
unique() %>%
dplyr::count(Kids_First_Participant_ID) %>%
# dplyr::mutate(kids_dec_n = glue::glue("{Kids_First_Participant_ID} (N={n})")) %>%
dplyr::rename(kids_deceased_bs_number = n) %>%
filter(!kids_deceased_bs_number <= 1) %>%
write_tsv(file.path(results_dir, "kids_dec_multiple_bs_list.tsv"))
###-----------------------------------------------------
# let's look into PT_3KM9W8S8
PT_3KM9W8S8_df <- nautilus_dec_df %>%
filter(Kids_First_Participant_ID == "PT_3KM9W8S8",
Kids_First_Biospecimen_ID == "BS_2NQXY528") %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, `Note field from Nautilus of initial parent`) %>%
unique() %>%
write_tsv(file.path(results_dir, "nautilus_dec_list.tsv"))
#bs_id <- print(unique(PT_3KM9W8S8_df$Kids_First_Biospecimen_ID))
Create pyclone input files
We need to generate the input files according to the method’s
template. Phylogenetic methods require at least 2 samples per tumor site
(multiregional sampling per anatomical site). Here, we will consider
kids samples with more than 2 timepoints with one or more biospecimens.
We will compare later differences in samples with single vs multiple
biospecimens.
# Create pyclone df for all samples
pyclone_all_samples_df <- nautilus_dec_df %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, cg_id,
cg_id_kids, mutation_id, sample_id,
tumor_descriptor, Chromosome, Start_Position,
Reference_Allele, Tumor_Seq_Allele2, t_ref_count, t_alt_count,
normal_cn, tumor_fraction, mutation_count) %>%
# change names to match input requirements
dplyr::rename("ref_counts" = "t_ref_count",
"alt_counts" = "t_alt_count",
"tumour_content" = "tumor_fraction") %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, tumor_descriptor,
cg_id_kids, mutation_id, sample_id, ref_counts,
alt_counts, normal_cn, tumour_content, mutation_count)
# I will test one HGG dataset for now
# PT_Z4BF2NSB
#PT_Z4BF2NSB_DF <- pyclone_all_samples_df %>%
# filter(Kids_First_Participant_ID == "PT_Z4BF2NSB")
Add major_cn, minor_cn columns from cns files
# Create and save pyclone_input!
pyclone_input <- data_list_bind %>%
# Rename to match input format
# To figure out if the assignment is correct
dplyr::rename("major_cn" = "cn1",
"minor_cn" = "cn2") %>%
left_join(pyclone_all_samples_df) %>%
filter(!is.na(major_cn),
!is.na(minor_cn)) %>%
dplyr::mutate(mutation_id = paste(mutation_id, major_cn, minor_cn,sep = ":")) %>%
select(Kids_First_Participant_ID, tumor_descriptor, cg_id_kids, mutation_id, sample_id, ref_counts,
alt_counts, normal_cn, major_cn, minor_cn, tumour_content, mutation_count) %>%
# To ensure there are no duplicated entries in the dataframe
distinct()
Joining with `by = join_by(Kids_First_Biospecimen_ID)`Warning: Detected an unexpected many-to-many relationship between `x` and `y`.
Plot depth coverage
# Read color palette
palette_df <- readr::read_tsv(palette_file, guess_max = 100000, show_col_types = FALSE) %>%
mutate(tumor_descriptor = color_names)
# Define and order palette
palette <- palette_df$hex_codes
names(palette) <- palette_df$tumor_descriptor
# Define timepoints
timepoints = c("Diagnosis", "Progressive", "Recurrence", "Deceased", "Second Malignancy", "Unavailable")
for (i in 1:length(data_dir) ) {
# Create sample_name
sample_name <- unique(as.character(gsub(".call.cns", "", str_split_fixed(data_dir[i], "/", 13)[,11])))
sample_name <- sort(sample_name, decreasing = FALSE)
print(sample_name)
for (x in seq_along(sample_name) ) {
pyclone_input_subset <- pyclone_input %>%
filter(Kids_First_Participant_ID == sample_name[x]) %>%
select(-c(Kids_First_Participant_ID)) #%>%
# mutate(tumor_descriptor = factor(tumor_descriptor),
# tumor_descriptor = fct_relevel(tumor_descriptor, timepoints)) %>%
#arrange(tumor_descriptor, sample_id)
# Make this reproducible
set.seed(2023)
# Define label for plots
Timepoint <- factor(x = pyclone_input_subset$tumor_descriptor, levels = timepoints)
#######################
# Create bxp ref_counts
p <- print(ggplot(pyclone_input_subset, aes(sample_id, ref_counts, color = Timepoint)) +
geom_jitter(width = 0.15, size = 0.7, alpha = 0.6) +
ggplot2::geom_boxplot(color = "black",
size = 0.25,
alpha = 0,
coef = 0) + # remove whiskers
theme_Publication() +
scale_color_manual(values = palette,
breaks = sort(names(palette))) +
#rotate() +
theme(axis.text.x = element_text(angle = 90)) +
stat_summary(fun.y=mean,shape=1,col='black',geom='point') +
labs(title = sample_name[x],
x = "sample_id",
y = "ref_counts",
color = "Timepoint"))
# Save the plot
ggsave(filename = paste0(sample_name[x], "-ref_counts.pdf"),
path = pyclone_plots_dir,
width = 6,
height = 5,
device = "pdf",
useDingbats = FALSE)
#######################
# Create bxp alt_counts
p <- print(ggplot(pyclone_input_subset, aes(sample_id, alt_counts, color = Timepoint)) +
geom_jitter(width = 0.15, size = 0.7, alpha = 0.6) +
ggplot2::geom_boxplot(color = "black",
size = 0.25,
alpha = 0,
coef = 0) + # remove whiskers
theme_Publication() +
scale_color_manual(values = palette,
breaks = sort(names(palette))) +
#rotate() +
theme(axis.text.x = element_text(angle = 90)) +
stat_summary(fun.y=mean,shape=1,col='black',geom='point') +
labs(title = sample_name[x],
x = "sample_id",
y = "alt_counts",
color = "Timepoint"))
# Save the plot
ggsave(filename = paste0(sample_name[x], "-alt_counts.pdf"),
path = pyclone_plots_dir,
width = 6,
height = 5,
device = "pdf",
useDingbats = FALSE)
}
}
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
























































Total number of mutations across timepoints and biospecimen sample
per Patient case
for (i in 1:length(data_dir) ) {
# Create sample_name
sample_name <- unique(as.character(gsub(".call.cns", "", str_split_fixed(data_dir[i], "/", 13)[,11])))
sample_name <- sort(sample_name, decreasing = FALSE)
print(sample_name)
for (x in seq_along(sample_name) ) {
pyclone_input_subset <- pyclone_input %>%
filter(Kids_First_Participant_ID == sample_name[x]) %>%
select(-c(Kids_First_Participant_ID))
# Define parameters for function
ylim = max(pyclone_input_subset$mutation_count)
# Rename legend for timepoints
Timepoint <- factor(pyclone_input_subset$tumor_descriptor)
# Plot stacked barplot
print(ggplot(pyclone_input_subset, aes(x = sample_id,
y = mutation_count,
fill = Timepoint)) +
geom_col(position = position_stack(reverse = TRUE)) +
geom_bar(stat = "identity", width = 0.5) +
scale_fill_manual(values = palette, breaks = sort(names(palette))) +
theme_Publication() +
theme(axis.text.x = element_text(angle = 85,
hjust = 1,
vjust = 1)) +
labs(title = paste(sample_name)) +
labs(x = "sample_id", y = "Total Mutations") +
ylim(0, ylim))
}
}
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"




























sessionInfo()
R version 4.3.1 (2023-06-16)
Platform: x86_64-apple-darwin20 (64-bit)
Running under: macOS Ventura 13.6.1
Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/lib/libRlapack.dylib; LAPACK version 3.11.0
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
time zone: America/New_York
tzcode source: internal
attached base packages:
[1] grid stats graphics grDevices utils datasets methods base
other attached packages:
[1] ggthemes_4.2.4 cDriver_0.4.2 readxl_1.4.3 gridExtra_2.3 clonevol_0.99.11 lubridate_1.9.3 forcats_1.0.0 stringr_1.5.1 dplyr_1.1.4
[10] purrr_1.0.2 readr_2.1.4 tidyr_1.3.0 tibble_3.2.1 ggplot2_3.4.4 tidyverse_2.0.0 dndscv_0.0.1.0
loaded via a namespace (and not attached):
[1] ade4_1.7-22 tidyselect_1.2.0 Rmpfr_0.9-3 farver_2.1.1 Biostrings_2.70.1 bitops_1.0-7
[7] fastmap_1.1.1 RCurl_1.98-1.13 digest_0.6.33 timechange_0.2.0 lifecycle_1.0.4 magrittr_2.0.3
[13] compiler_4.3.1 rlang_1.1.2 sass_0.4.7 tools_4.3.1 utf8_1.2.4 yaml_2.3.7
[19] data.table_1.14.8 knitr_1.45 labeling_0.4.3 bit_4.0.5 abind_1.4-5 BiocParallel_1.36.0
[25] withr_2.5.2 BiocGenerics_0.48.1 stats4_4.3.1 fansi_1.0.5 colorspace_2.1-0 scales_1.2.1
[31] MASS_7.3-60 cli_3.6.1 rmarkdown_2.25 crayon_1.5.2 ragg_1.2.6 generics_0.1.3
[37] rstudioapi_0.15.0 tzdb_0.4.0 cachem_1.0.8 zlibbioc_1.48.0 parallel_4.3.1 cellranger_1.1.0
[43] XVector_0.42.0 vctrs_0.6.4 jsonlite_1.8.7 carData_3.0-5 car_3.1-2 IRanges_2.36.0
[49] hms_1.1.3 S4Vectors_0.40.1 bit64_4.0.5 rstatix_0.7.2 seqinr_4.2-30 systemfonts_1.0.5
[55] jquerylib_0.1.4 glue_1.6.2 codetools_0.2-19 stringi_1.8.1 gtable_0.3.4 GenomeInfoDb_1.38.1
[61] GenomicRanges_1.54.1 gmp_0.7-2 munsell_0.5.0 pillar_1.9.0 htmltools_0.5.7 GenomeInfoDbData_1.2.11
[67] R6_2.5.1 textshaping_0.3.7 rprojroot_2.0.4 vroom_1.6.4 evaluate_0.23 backports_1.4.1
[73] Rsamtools_2.18.0 broom_1.0.5 bslib_0.5.1 Rcpp_1.0.11 xfun_0.41 pkgconfig_2.0.3
LS0tCnRpdGxlOiAiSW5mZXJlbmNlIG9mIHN1YmNsb25hbCBhcmNoaXRlY3R1cmUgb2YgdHVtb3JzIGFjcm9zcyBtdWx0aXBsZSB0aW1lcG9pbnRzIGluIHRoZSBwYWlyZWQgbG9uZ2l0dWRpbmFsIChQTCkgY29ob3J0IgphdXRob3I6ICdBbnRvbmlhIENocm9uaSA8Y2hyb25pYUBjaG9wLmVkdT4gZm9yIEQzQicKZGF0ZTogIjIwMjMiCm91dHB1dDoKICBodG1sX25vdGVib29rOgogICAgdG9jOiBUUlVFCiAgICB0b2NfZmxvYXQ6IFRSVUUKLS0tCgojIEJhY2tncm91bmQKClRoaXMgbm90ZWJvb2sgaXMgZm9yIHByZXBhcmluZyBpbnB1dCBmaWxlcyBmb3IgW3B5Y2xvbmUtdmldKGh0dHBzOi8vZ2l0aHViLmNvbS9Sb3RoLUxhYi9weWNsb25lLXZpKS4gVGhlIHNhbWUgZmlsZXMgd2l0aCBhIG1pbm9yIG1vZGlmaWNhdGlvbiB3aWxsIGJlIHVzZWQgYXMgYW4gaW5wdXQgZm9yIFtGYXN0Q2xvbmVdKGh0dHBzOi8vZ2l0aHViLmNvbS9HdWFuTGFiL0Zhc3RDbG9uZV9HdWFuTGFiKS4KCgojIFNldCB1cAoKYGBge3IgbG9hZC1saWJyYXJ5fQpzdXBwcmVzc1BhY2thZ2VTdGFydHVwTWVzc2FnZXMoewogIGxpYnJhcnkodGlkeXZlcnNlKQogIGxpYnJhcnkocmVhZHhsKQogIGxpYnJhcnkoY0RyaXZlcikgIyBDYWxjdWxhdGUgQ0NGLCBodHRwczovL2dpdGh1Yi5jb20vaGFuYXN1c2FrL2NEcml2ZXIvCn0pCmBgYAoKIyBEaXJlY3RvcmllcyBhbmQgRmlsZSBJbnB1dHMvT3V0cHV0cwoKYGBge3Igc2V0LWRpci1hbmQtZmlsZS1uYW1lc30KIyBEZXRlY3QgdGhlICIuZ2l0IiBmb2xkZXIgLS0gdGhpcyB3aWxsIGJlIGluIHRoZSBwcm9qZWN0IHJvb3QgZGlyZWN0b3J5CiMgVXNlIHRoaXMgYXMgdGhlIHJvb3QgZGlyZWN0b3J5IHRvIGVuc3VyZSBwcm9wZXIgc291cmNpbmcgb2YgZnVuY3Rpb25zIAojIG5vIG1hdHRlciB3aGVyZSB0aGlzIGlzIGNhbGxlZCBmcm9tCnJvb3RfZGlyIDwtIHJwcm9qcm9vdDo6ZmluZF9yb290KHJwcm9qcm9vdDo6aGFzX2RpcigiLmdpdCIpKQphbmFseXNpc19kaXIgPC0gZmlsZS5wYXRoKHJvb3RfZGlyLCAiYW5hbHlzZXMiLCAidHVtb3ItY2xvbmUtaW5mZXJlbmNlIikKaW5wdXRfZGlyIDwtIGZpbGUucGF0aChhbmFseXNpc19kaXIsICJpbnB1dCIpCmRhdGFfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImRhdGEiKQpmaWxlc19kaXIgPC0gZmlsZS5wYXRoKHJvb3RfZGlyLCAiYW5hbHlzZXMiLCAic2FtcGxlLWRpc3RyaWJ1dGlvbi1hbmFseXNpcyIsICJyZXN1bHRzIikKbWFmX2ZpbGVzX2RpciA8LSBmaWxlLnBhdGgocm9vdF9kaXIsICJhbmFseXNlcyIsICJ0bWItdmFmLWxvbmdpdHVkaW5hbCIsICJyZXN1bHRzIikKCiMgSW5wdXQgZmlsZXMKI3BidGFfZmlsZSA8LSBmaWxlLnBhdGgoZmlsZXNfZGlyLCAicGJ0YS50c3YiKSAjIGZpbGUgZnJvbSBhZGQtc2FtcGxlLWRpc3RyaWJ1dGlvbiBtb2R1bGUKbWFmX2ZpbGUgPC0gZmlsZS5wYXRoKG1hZl9maWxlc19kaXIsICJ0bWJfdmFmX2dlbm9taWMudHN2IikKZ2Vub21pY19wYWlyZWRfZmlsZSA8LSBmaWxlLnBhdGgoZmlsZXNfZGlyLCAiZ2Vub21pY19hc3NheXNfbWF0Y2hlZF90aW1lX3BvaW50cy50c3YiKSAjIGZpbGUgZnJvbSBhZGQtc2FtcGxlLWRpc3RyaWJ1dGlvbiBtb2R1bGUKbmF1dGlsdXNfZGVjX2ZpbGUgPC0gZmlsZS5wYXRoKGlucHV0X2RpciwgImRlY2Vhc2VkX3NhbXBsZXMueGxzeCIpIApwYWxldHRlX2ZpbGUgPC0gZmlsZS5wYXRoKHJvb3RfZGlyLCAiZmlndXJlcyIsICJwYWxldHRlcyIsICJ0dW1vcl9kZXNjcmlwdG9yX2NvbG9yX3BhbGV0dGUudHN2IikKCgojIEVYQU1QTEUgCiMgQnV0IHRvIHJlcGxhY2Ugd2l0aCBkaXIgd2l0aCBjbnMgZGF0YSBpbmZlcnJlZCBieSBDTlZraXQKIyBQVF9aNEJGMk5TQl9kaXIgPC0gZmlsZS5wYXRoKGlucHV0X2RpciwgImNudmtpdF9kYXRhX2V4YW1wbGUvUFRfWjRCRjJOU0IiKSAKY25zX2RpciA8LSBmaWxlLnBhdGgoaW5wdXRfZGlyLCAiY252a2l0X2RhdGEiKSAKCiMgRmlsZSBwYXRoIHRvIHJlc3VsdHMgZGlyZWN0b3J5CnJlc3VsdHNfZGlyIDwtCiAgZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgInJlc3VsdHMiKQppZiAoIWRpci5leGlzdHMocmVzdWx0c19kaXIpKSB7CiAgZGlyLmNyZWF0ZShyZXN1bHRzX2RpcikKfQoKIyBGaWxlIHBhdGggdG8gaW5wdXQgZGlyZWN0b3J5CnB5Y2xvbmV2aV9pbnB1dF9kaXIgPC0KICBmaWxlLnBhdGgoYW5hbHlzaXNfZGlyLCAicmVzdWx0cyIsICJweWNsb25lLXZpLWlucHV0IikKaWYgKCFkaXIuZXhpc3RzKHB5Y2xvbmV2aV9pbnB1dF9kaXIpKSB7CiAgZGlyLmNyZWF0ZShweWNsb25ldmlfaW5wdXRfZGlyKQp9CgojIEZpbGUgcGF0aCB0byBpbnB1dCBkaXJlY3RvcnkKZmFzdGNsb25lX2lucHV0X2RpciA8LQogIGZpbGUucGF0aChhbmFseXNpc19kaXIsICJyZXN1bHRzIiwgImZhc3RjbG9uZS1pbnB1dCIpCmlmICghZGlyLmV4aXN0cyhweWNsb25ldmlfaW5wdXRfZGlyKSkgewogIGRpci5jcmVhdGUocHljbG9uZXZpX2lucHV0X2RpcikKfQoKCiMgRmlsZSBwYXRoIHRvIHBsb3QgZGlyZWN0b3J5CnB5Y2xvbmVfcGxvdHNfZGlyIDwtCiAgZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgInBsb3RzIiwgInB5Y2xvbmUtdmkiKQppZiAoIWRpci5leGlzdHMocHljbG9uZV9wbG90c19kaXIpKSB7CiAgZGlyLmNyZWF0ZShweWNsb25lX3Bsb3RzX2RpcikKfQoKCnNvdXJjZShwYXN0ZTAocm9vdF9kaXIsICIvZmlndXJlcy9zY3JpcHRzL3RoZW1lLlIiKSkKYGBgCgojIExvYWQgYW5kIHByb2Nlc3MgZGF0YQoKYGBge3IgbG9hZC1wcm9jZXNzLWlucHV0c30KIyBSZWFkIG1hZgptYWZfZGYgPC0gcmVhZHI6OnJlYWRfdHN2KG1hZl9maWxlLCBndWVzc19tYXggPSAxMDAwMDAsIHNob3dfY29sX3R5cGVzID0gRkFMU0UpICU+JQogIHNlbGVjdChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lELCBjZ19pZCwgdHVtb3JfZGVzY3JpcHRvciwKICAgICAgICAgc2FtcGxlX2lkLCBhbGlxdW90X2lkLCBDaHJvbW9zb21lLCAKICAgICAgICAgU3RhcnRfUG9zaXRpb24sIFJlZmVyZW5jZV9BbGxlbGUsIFR1bW9yX1NlcV9BbGxlbGUxLCBUdW1vcl9TZXFfQWxsZWxlMiwKICAgICAgICAgdF9yZWZfY291bnQsIHRfYWx0X2NvdW50LCB0bWIsIAogICAgICAgICB0dW1vcl9mcmFjdGlvbiwgdHVtb3JfcGxvaWR5LCBWQUYsIEh1Z29fU3ltYm9sLAogICAgICAgICBtdXRhdGlvbl9jb3VudCkgJT4lIAoKICAjIFJlbW92ZSBoeXBlcm11dGFudHMKICBmaWx0ZXIoIXRtYiA+PSAxMCwgCiAgCiAgIyBUaGVyZSBhcmUgYWx0ZXJhdGlvbnMgd2l0aCBoaWdoIG51bWJlciBvZiByZWFkIGNvdW50cy4KICAjIFdlIHdpbGwgZXhjbHVkZSB0aG9zZSB3aXRoID4xMDAwIGZvciBub3cuCiAgICAgICAgICF0X2FsdF9jb3VudCA+PSAxMDAwLAogICAgICAgICAhdF9yZWZfY291bnQgPj0gMTAwMCkgICU+JSAKICAKICAjIEFkZCBgbm9ybWFsX2NuYDogVG90YWwgY29weSBudW1iZXIgb2Ygc2VnbWVudCBpbiBoZWFsdGh5IHRpc3N1ZS4gCiAgIyBGb3IgYXV0b3NvbWUgdGhpcyB3aWxsIGJlIHR3byBhbmQgbWFsZSBzZXggY2hyb21vc29tZXMgb25lLgogICMgU2VlOiBodHRwczovL2dpdGh1Yi5jb20vUm90aC1MYWIvcHljbG9uZS12aQogIG11dGF0ZShub3JtYWxfY24gPSBjYXNlX3doZW4oZ3JlcGwoImNoclkiLCBDaHJvbW9zb21lKSB+ICIxIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBUUlVFIH4gIjIiKSkKIyBDYWxjdWxhdGUgQ0NGCm1hZl9kZiA8LSBDQ0YobWFmX2RmKQogCiMgUmVhZCBwYXRpZW50IGxpc3QKZ2Vub21pY19wYWlyZWRfZGYgPC0gcmVhZHI6OnJlYWRfdHN2KGdlbm9taWNfcGFpcmVkX2ZpbGUsIGd1ZXNzX21heCA9IDEwMDAwMCwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkgJT4lIAogIHNlbGVjdCghYyhjYW5jZXJfZ3JvdXAsIGV4cGVyaW1lbnRhbF9zdHJhdGVneSkpICU+JSAKICBsZWZ0X2pvaW4obWFmX2RmLCBieSA9IGMoIktpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQiKSkgJT4lICMgY3JlYXRlIHVuaXF1ZSBpZGVudGlmaWVycwogIGRwbHlyOjptdXRhdGUobXV0YXRpb25faWQgPSBwYXN0ZShLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCB0dW1vcl9kZXNjcmlwdG9yLCBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lELCBDaHJvbW9zb21lLCBTdGFydF9Qb3NpdGlvbiwgUmVmZXJlbmNlX0FsbGVsZSwgVHVtb3JfU2VxX0FsbGVsZTIsIHNlcCA9ICI6IiksCiAgICAgICAgICAgICAgICBjZ19pZF9raWRzID0gcGFzdGUoY2dfaWQsIEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIHNlcCA9ICJfIiksCiAgICAgICAgICAgICAgICBjZ19pZF9raWRzID0gc3RyX3JlcGxhY2UoY2dfaWRfa2lkcywgIi98LSIsICJfIiksCiAgICAgICAgICAgICAgICBjZ19pZF9raWRzID0gc3RyX3JlcGxhY2VfYWxsKGNnX2lkX2tpZHMsICIgIiwgIl8iKSwKICAgICAgICAgICAgICAgIGNnX2lkID0gc3RyX3JlcGxhY2UoY2dfaWQsICIvfC0iLCAiXyIpLAogICAgICAgICAgICAgICAgY2dfaWQgPSBzdHJfcmVwbGFjZV9hbGwoY2dfaWQsICIgIiwgIl8iKSwKICAgICAgICAgICAgICAgIHNhbXBsZV9pZCA9IHBhc3RlKHR1bW9yX2Rlc2NyaXB0b3IsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIHNlcCA9ICI6IikpCgoKIyBMZXQncyBjb3VudCAjc3BlY2ltZW5zIHBlciBzYW1wbGUgCiMgV2Ugd2lsbCByZW1vdmUgYW55IHNhbXBsZXMgd2l0aCBsZXNzIHRoYW4gMiBzcGVjaW1lbnMgYXMgcGh5bG9nZW5pZXMgcmVxdWlyZSBhdCBsZWFzdCAzIHRheGEKa2lkc19zcGVjaW1lbnNfbl9kZiA8LSBnZW5vbWljX3BhaXJlZF9kZiAlPiUgCiAgc2VsZWN0KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIHR1bW9yX2Rlc2NyaXB0b3IpICU+JSAKICB1bmlxdWUoKSAlPiUgCiAgZHBseXI6OmNvdW50KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpICU+JSAKICBkcGx5cjo6bXV0YXRlKGtpZHNfc3BlY2ltZW5zX24gPSBnbHVlOjpnbHVlKCJ7S2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRH0gIChOPXtufSkiKSkgJT4lCiAgZHBseXI6OnJlbmFtZShraWRzX3NwZWNpbWVuc19udW1iZXIgPSBuKSAlPiUgCiAgZmlsdGVyKCFraWRzX3NwZWNpbWVuc19udW1iZXIgPD0gMikgJT4lCgogIGxlZnRfam9pbihnZW5vbWljX3BhaXJlZF9kZiwgYnkgPSBjKCJLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEIikpICU+JQogICNsZWZ0X2pvaW4obWFmX2RmKSAlPiUgCiAgZmlsdGVyKCFpcy5uYShDaHJvbW9zb21lKSkgCiAgCiMgTGV0J3MgY29uZmlybSBhbmQgYWRkIHRoZSBudW1iZXIgb2YgdGltZXBvaW50cyBwZXIgc2FtcGxlCiMgSW4gY2FzZSB3ZSBsb3N0IGFueSBmcm9tIGZpbHRlcmluZyBoeXBlcm11dGFudHMgYW5kIGhpZ2ggcmVhZHMvYWx0ZXJhdGlvbgp0aW1lcG9pbnRzX25fZGYgPC0ga2lkc19zcGVjaW1lbnNfbl9kZiAlPiUgCiAgc2VsZWN0KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIHR1bW9yX2Rlc2NyaXB0b3IpICU+JSAKICB1bmlxdWUoKSAlPiUgCiAgZHBseXI6OmNvdW50KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpICU+JSAKICBkcGx5cjo6bXV0YXRlKGtpZHNfdGltZXBvaW50c19uID0gZ2x1ZTo6Z2x1ZSgie0tpZHNfRmlyc3RfUGFydGljaXBhbnRfSUR9ICAoTj17bn0pIikpICU+JQogIGRwbHlyOjpyZW5hbWUoa2lkc190aW1lcG9pbnRzX251bWJlciA9IG4pICU+JSAKICBmaWx0ZXIoIWtpZHNfdGltZXBvaW50c19udW1iZXIgPD0gMSkKCmRmIDwtIHRpbWVwb2ludHNfbl9kZiAlPiUgCiAgbGVmdF9qb2luKGtpZHNfc3BlY2ltZW5zX25fZGYsIGJ5ID0gYygiS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCIpKSAlPiUgCiAgd3JpdGVfdHN2KGZpbGUucGF0aChyZXN1bHRzX2RpciwgInNhbXBsZXNfZWxpZ2libGVfZm9yX3BoeWxvZ2VueS50c3YiKSkKCiMgTGlzdCB3aXRoIHNhbXBsZXMgZWxpZ2libGUgZm9yIHBoeWxvZ2VueQojIEkgYWRkZWQgdGhlIGluZm9ybWF0aW9uIGFib3V0IHRvIHVzZSBmb3IgcGh5bG9nZW5ldGljIGluZmVyZW5jZXMKbGlzdF9kZiA8LSBkZiAlPiUgCiAgc2VsZWN0KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpICU+JSAKICB1bmlxdWUoKSAlPiUgCiAgbXV0YXRlKHNvbWF0aWNfZ2VybWxpbmVfcGh5bG9nZW55ID0gY2FzZV93aGVuKGdyZXBsKCJQVF9LWjU2WEhKVHxQVF9LVFJKOFRGWSIsIEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpIH4gInllcyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBUUlVFIH4gIm5vdF95ZXQiKSkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJzYW1wbGVzX2VsaWdpYmxlX2Zvcl9waHlsb2dlbnlfbGlzdC50c3YiKSkKCmBgYAoKIyMgQWRkIE5hdXRpbHVzIGxvY2F0aW9uIGZvciBEZWNlYXNlZCBzcGVjaW1lbnMKCmBgYHtyIGFkZC1OYXV0aWx1cy1sb2NhdGlvbn0KbmF1dGlsdXNfZGVjX2RmIDwtIHJlYWRfZXhjZWwobmF1dGlsdXNfZGVjX2ZpbGUpICU+JSAKICByaWdodF9qb2luKGRmLCBieSA9IGMoInNhbXBsZV9pZCIsICJhbGlxdW90X2lkIiwgInR1bW9yX2Rlc2NyaXB0b3IiKSkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJuYXV0aWx1c19kZWMudHN2IikpCgpsaXN0X2RmIDwtIG5hdXRpbHVzX2RlY19kZiAlPiUgCiAgc2VsZWN0KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIGBOb3RlIGZpZWxkIGZyb20gTmF1dGlsdXMgb2YgaW5pdGlhbCBwYXJlbnRgKSAlPiUgCiAgdW5pcXVlKCkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJuYXV0aWx1c19kZWNfbGlzdC50c3YiKSkKCmBgYAoKIyMgT3RoZXIsIG1heWJlIHRvIGRlbGV0ZQoKYGBgIHtyIG90aGVyLXRoaW5nc30KCiMgTWFrZSBsaXN0IHdpdGggc2FtcGxlcyB3aXRoID4yIGJpb3NwZWNpbWVucyBhdCBEZWNlYXNlZCB0aW1lcG9pbnQKZGVjX25fZGYgPC0gZGYgJT4lIAogIGZpbHRlcih0dW1vcl9kZXNjcmlwdG9yID09ICJEZWNlYXNlZCIpICU+JSAKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgdHVtb3JfZGVzY3JpcHRvciwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCkgJT4lIAogIHVuaXF1ZSgpICU+JSAKICBkcGx5cjo6Y291bnQoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkgJT4lIAogIyBkcGx5cjo6bXV0YXRlKGtpZHNfZGVjX24gPSBnbHVlOjpnbHVlKCJ7S2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRH0gIChOPXtufSkiKSkgJT4lCiAgZHBseXI6OnJlbmFtZShraWRzX2RlY2Vhc2VkX2JzX251bWJlciA9IG4pICU+JSAKICBmaWx0ZXIoIWtpZHNfZGVjZWFzZWRfYnNfbnVtYmVyIDw9IDEpICU+JSAKICB3cml0ZV90c3YoZmlsZS5wYXRoKHJlc3VsdHNfZGlyLCAia2lkc19kZWNfbXVsdGlwbGVfYnNfbGlzdC50c3YiKSkKCiMjIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCiMgbGV0J3MgbG9vayBpbnRvIFBUXzNLTTlXOFM4ClBUXzNLTTlXOFM4X2RmIDwtIG5hdXRpbHVzX2RlY19kZiAlPiUgCiAgZmlsdGVyKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQgPT0gIlBUXzNLTTlXOFM4IiwKICAgICAgICAgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCA9PSAiQlNfMk5RWFk1MjgiKSAlPiUgCiAgc2VsZWN0KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIGBOb3RlIGZpZWxkIGZyb20gTmF1dGlsdXMgb2YgaW5pdGlhbCBwYXJlbnRgKSAlPiUgCiAgdW5pcXVlKCkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJuYXV0aWx1c19kZWNfbGlzdC50c3YiKSkKICAKI2JzX2lkIDwtIHByaW50KHVuaXF1ZShQVF8zS005VzhTOF9kZiRLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEKSkKYGBgCgoKIyBDcmVhdGUgcHljbG9uZSBpbnB1dCBmaWxlcyAKCldlIG5lZWQgdG8gZ2VuZXJhdGUgdGhlIGlucHV0IGZpbGVzIGFjY29yZGluZyB0byB0aGUgbWV0aG9kJ3MgdGVtcGxhdGUuIApQaHlsb2dlbmV0aWMgbWV0aG9kcyByZXF1aXJlIGF0IGxlYXN0IDIgc2FtcGxlcyBwZXIgdHVtb3Igc2l0ZSAobXVsdGlyZWdpb25hbCBzYW1wbGluZyBwZXIgYW5hdG9taWNhbCBzaXRlKS4KSGVyZSwgd2Ugd2lsbCBjb25zaWRlciBraWRzIHNhbXBsZXMgd2l0aCBtb3JlIHRoYW4gMiB0aW1lcG9pbnRzIHdpdGggb25lIG9yIG1vcmUgYmlvc3BlY2ltZW5zLiAKV2Ugd2lsbCBjb21wYXJlIGxhdGVyIGRpZmZlcmVuY2VzIGluIHNhbXBsZXMgd2l0aCBzaW5nbGUgdnMgbXVsdGlwbGUgYmlvc3BlY2ltZW5zLgoKYGBge3IgY3JlYXRlLXB5Y2xvbmUtYWxsLXNhbXBsZXMtZGZ9CiMgQ3JlYXRlIHB5Y2xvbmUgZGYgZm9yIGFsbCBzYW1wbGVzCnB5Y2xvbmVfYWxsX3NhbXBsZXNfZGYgPC0gbmF1dGlsdXNfZGVjX2RmICU+JQogIHNlbGVjdChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lELCBjZ19pZCwgCiAgICAgICAgIGNnX2lkX2tpZHMsIG11dGF0aW9uX2lkLCBzYW1wbGVfaWQsCiAgICAgICAgIHR1bW9yX2Rlc2NyaXB0b3IsIENocm9tb3NvbWUsIFN0YXJ0X1Bvc2l0aW9uLCAKICAgICAgICAgUmVmZXJlbmNlX0FsbGVsZSwgVHVtb3JfU2VxX0FsbGVsZTIsIHRfcmVmX2NvdW50LCB0X2FsdF9jb3VudCwKICAgICAgICAgbm9ybWFsX2NuLCB0dW1vcl9mcmFjdGlvbiwgbXV0YXRpb25fY291bnQpICU+JSAKICAKICAjIGNoYW5nZSBuYW1lcyB0byBtYXRjaCBpbnB1dCByZXF1aXJlbWVudHMKICBkcGx5cjo6cmVuYW1lKCJyZWZfY291bnRzIiA9ICJ0X3JlZl9jb3VudCIsIAogICAgICAgICAgICAgICAgImFsdF9jb3VudHMiID0gInRfYWx0X2NvdW50IiwKICAgICAgICAgICAgICAgICJ0dW1vdXJfY29udGVudCIgPSAidHVtb3JfZnJhY3Rpb24iKSAlPiUKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgdHVtb3JfZGVzY3JpcHRvciwgCiAgICAgICAgIGNnX2lkX2tpZHMsIG11dGF0aW9uX2lkLCBzYW1wbGVfaWQsIHJlZl9jb3VudHMsIAogICAgICAgICBhbHRfY291bnRzLCBub3JtYWxfY24sIHR1bW91cl9jb250ZW50LCBtdXRhdGlvbl9jb3VudCkKIAoKIyBJIHdpbGwgdGVzdCBvbmUgSEdHIGRhdGFzZXQgZm9yIG5vdwojIFBUX1o0QkYyTlNCCgojUFRfWjRCRjJOU0JfREYgPC0gcHljbG9uZV9hbGxfc2FtcGxlc19kZiAlPiUgCiMgIGZpbHRlcihLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEID09ICJQVF9aNEJGMk5TQiIpCiAgICAgICAgICAgICAgICAgICAKYGBgCgojIyBBZGQgbWFqb3JfY24sIG1pbm9yX2NuIGNvbHVtbnMgZnJvbSBjbnMgZmlsZXMKCmBgYHtyIHByb2Nlc3MtY25zLWZpbGVzfQpkYXRhX2RpciA8LSBkaXIocGF0aCA9IGNuc19kaXIsICBwYXR0ZXJuID0gIi5jYWxsLmNucyIsIGZ1bGwubmFtZXMgPSBUUlVFLCByZWN1cnNpdmUgPSBUUlVFKQoKIyBDcmVhdGUgbGlzdCAKZGF0YV9saXN0IDwtIGxpc3QoKQoKZm9yIChpIGluIDE6bGVuZ3RoKGRhdGFfZGlyKSApIHsgCiAgCiAgIyBDcmVhdGUgc2FtcGxlX25hbWUKICBzYW1wbGVfbmFtZSA8LSAgdW5pcXVlKGFzLmNoYXJhY3Rlcihnc3ViKCIuY2FsbC5jbnMiLCAiIiwgc3RyX3NwbGl0X2ZpeGVkKGRhdGFfZGlyW2ldLCAiLyIsIDEzKVssMTFdKSkpCiAgc2FtcGxlX25hbWUgPC0gc29ydChzYW1wbGVfbmFtZSwgZGVjcmVhc2luZyA9IEZBTFNFKQogIHByaW50KHNhbXBsZV9uYW1lKQogIAogIGZvciAoeCBpbiBzZXFfYWxvbmcoc2FtcGxlX25hbWUpICkgeyAKICAgIAogICAgZGF0YV9saXN0W1tpXV0gPC0gcmVhZC5jc3YoZGF0YV9kaXJbaV0sIGhlYWRlcj1ULCBzZXA9Ilx0IikgCiAgCiAgICAKICAjIENyZWF0ZSBmaWxlX25hbWUKICBmaWxlX25hbWUgPC0gZ3N1YigiLmNhbGwuY25zIiwgIiIsIHN0cl9zcGxpdF9maXhlZChkYXRhX2RpcltpXSwgIi8iLCAxMylbLDEyXSkKICBwcmludChmaWxlX25hbWUpCiAgCiAgZGF0YV9saXN0W1tpXV0gPC0gZGF0YV9saXN0W1tpXV0gJT4lIAogICAgbXV0YXRlKEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQgPSBmaWxlX25hbWUpCiAgCiAgIyBUaGUgZm9sbG93aW5nIGNvZGUgYXNzaWducyBuYW1lIHRvIGRmCiAgIyBCdXQgYWRkcyBhbiBleHRyYSBkZiBpbnRvIHRoZSBsaXN0IC0gZG9udCBrbm93IHlldCB3aHkgYnV0IGxldCdzIHJlbW92ZSwgaXQncyBub3QgbmVjZXNzYXJ5CiAgI2RmIDwtIGFzc2lnbihmaWxlX25hbWUsIGRhdGFfbGlzdCkgCiAgI2RhdGFfbGlzdFtbZmlsZV9uYW1lXV0gPC0gZGYKICB9Cn0KCiMgQmluZCBhbGwgZGYgZnJvbSBsaXN0IApkYXRhX2xpc3RfYmluZCA8LSBkcGx5cjo6YmluZF9yb3dzKGRhdGFfbGlzdCkgCgojIENyZWF0ZSBhbmQgc2F2ZSBweWNsb25lX2lucHV0IQpweWNsb25lX2lucHV0IDwtIGRhdGFfbGlzdF9iaW5kICU+JSAKICAKICAjIFJlbmFtZSB0byBtYXRjaCBpbnB1dCBmb3JtYXQKICAjIFRvIGZpZ3VyZSBvdXQgaWYgdGhlIGFzc2lnbm1lbnQgaXMgY29ycmVjdAogIGRwbHlyOjpyZW5hbWUoIm1ham9yX2NuIiA9ICJjbjEiLCAKICAgICAgICAgICAgICAgICJtaW5vcl9jbiIgPSAiY24yIikgJT4lIAoKICBsZWZ0X2pvaW4ocHljbG9uZV9hbGxfc2FtcGxlc19kZikgJT4lIAogIGZpbHRlcighaXMubmEobWFqb3JfY24pLAogICAgICAgICAhaXMubmEobWlub3JfY24pKSAlPiUKICBkcGx5cjo6bXV0YXRlKG11dGF0aW9uX2lkID0gcGFzdGUobXV0YXRpb25faWQsIG1ham9yX2NuLCBtaW5vcl9jbixzZXAgPSAiOiIpKSAlPiUgCiAgc2VsZWN0KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIHR1bW9yX2Rlc2NyaXB0b3IsIGNnX2lkX2tpZHMsIG11dGF0aW9uX2lkLCBzYW1wbGVfaWQsIHJlZl9jb3VudHMsIAogICAgICAgICBhbHRfY291bnRzLCBub3JtYWxfY24sIG1ham9yX2NuLCBtaW5vcl9jbiwgdHVtb3VyX2NvbnRlbnQsIG11dGF0aW9uX2NvdW50KSAlPiUgCiAgCiAgIyBUbyBlbnN1cmUgdGhlcmUgYXJlIG5vIGR1cGxpY2F0ZWQgZW50cmllcyBpbiB0aGUgZGF0YWZyYW1lCiAgZGlzdGluY3QoKSAKCgpmb3IgKGkgaW4gMTpsZW5ndGgoZGF0YV9kaXIpICkgeyAKICAKICAjIENyZWF0ZSBzYW1wbGVfbmFtZQogIHNhbXBsZV9uYW1lIDwtICB1bmlxdWUoYXMuY2hhcmFjdGVyKGdzdWIoIi5jYWxsLmNucyIsICIiLCBzdHJfc3BsaXRfZml4ZWQoZGF0YV9kaXJbaV0sICIvIiwgMTMpWywxMV0pKSkKICBzYW1wbGVfbmFtZSA8LSBzb3J0KHNhbXBsZV9uYW1lLCBkZWNyZWFzaW5nID0gRkFMU0UpCiAgcHJpbnQoc2FtcGxlX25hbWUpCiAgCmZvciAoeCBpbiBzZXFfYWxvbmcoc2FtcGxlX25hbWUpICkgeyAKICAKICAjIFNhdmUgaW5wdXQgZmlsZSBmb3IgcHljbG9uZS12aQogIGZuYW1lIDwtIHBhc3RlMChweWNsb25ldmlfaW5wdXRfZGlyLCAiLyIsIHNhbXBsZV9uYW1lW3hdLCAiLnRzdiIpCiAgcHJpbnQoZm5hbWUpCiAgCiAgcHljbG9uZV9pbnB1dF9zdWJzZXQgPC0gcHljbG9uZV9pbnB1dCAlPiUKICAgIGZpbHRlcihLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEID09IHNhbXBsZV9uYW1lW3hdKSAlPiUgCiAgICAjc2VsZWN0KC1jKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpKSAlPiUgCiAgd3JpdGVfdHN2KGZpbGUucGF0aChmbmFtZSkpCiAgCiAgCiAgIyBTYXZlIGlucHV0IGZpbGUgZm9yIEZhc3RDbG9uZQogIGZhc3RjbG9uZV9mbmFtZSA8LSBwYXN0ZTAoZmFzdGNsb25lX2lucHV0X2RpciwgIi8iLCBzYW1wbGVfbmFtZVt4XSwgIi50c3YiKQogIHByaW50KGZhc3RjbG9uZV9mbmFtZSkKICAKICBmYXN0Y2xvbmVfaW5wdXRfc3Vic2V0IDwtIHB5Y2xvbmVfaW5wdXQgJT4lCiAgICBkcGx5cjo6cmVuYW1lKCJ2YXJfY291bnRzIiA9ICJhbHRfY291bnRzIikgJT4lIAogICAgZmlsdGVyKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQgPT0gc2FtcGxlX25hbWVbeF0pICU+JSAKICAgICNzZWxlY3QoLWMoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkpICU+JSAKICB3cml0ZV90c3YoZmlsZS5wYXRoKGZhc3RjbG9uZV9mbmFtZSkpCiAgICAKICAgIAogICAgCn0KfQoKCiMgVG8gZmluZCB0aGUgcG9zaXRpb24gb2YgZHVwbGljYXRlIGVsZW1lbnRzIGluIHgsIHVzZSB0aGlzOgojIGR1cGxpY2F0ZWQocHljbG9uZV9pbnB1dCkKCiMgRXh0cmFjdCBkdXBsaWNhdGUgZWxlbWVudHM6CiMgcHljbG9uZV9pbnB1dFtkdXBsaWNhdGVkKHB5Y2xvbmVfaW5wdXQpXQoKYGBgCgojIFBsb3QgZGVwdGggY292ZXJhZ2UKCmBgYHtyIGRlZmluZS1wYXJhbWV0ZXJzLWZvci1wbG90c30KIyBSZWFkIGNvbG9yIHBhbGV0dGUKcGFsZXR0ZV9kZiA8LSByZWFkcjo6cmVhZF90c3YocGFsZXR0ZV9maWxlLCBndWVzc19tYXggPSAxMDAwMDAsIHNob3dfY29sX3R5cGVzID0gRkFMU0UpICU+JSAKICBtdXRhdGUodHVtb3JfZGVzY3JpcHRvciA9IGNvbG9yX25hbWVzKQoKIyBEZWZpbmUgYW5kIG9yZGVyIHBhbGV0dGUKcGFsZXR0ZSA8LSBwYWxldHRlX2RmJGhleF9jb2RlcwpuYW1lcyhwYWxldHRlKSA8LSBwYWxldHRlX2RmJHR1bW9yX2Rlc2NyaXB0b3IKCiMgRGVmaW5lIHRpbWVwb2ludHMKdGltZXBvaW50cyA9IGMoIkRpYWdub3NpcyIsICJQcm9ncmVzc2l2ZSIsICJSZWN1cnJlbmNlIiwgIkRlY2Vhc2VkIiwgIlNlY29uZCBNYWxpZ25hbmN5IiwgIlVuYXZhaWxhYmxlIikKCmBgYAoKYGBgIHtyIGRlcHRoLWNvdmVyYWdlLCBmaWcud2lkdGggPSA2LCBmaWcuaGVpZ2h0ID0gNSwgZmlnLmZ1bGx3aWR0aCA9IFRSVUV9CmZvciAoaSBpbiAxOmxlbmd0aChkYXRhX2RpcikgKSB7IAogIAogICMgQ3JlYXRlIHNhbXBsZV9uYW1lCiAgc2FtcGxlX25hbWUgPC0gIHVuaXF1ZShhcy5jaGFyYWN0ZXIoZ3N1YigiLmNhbGwuY25zIiwgIiIsIHN0cl9zcGxpdF9maXhlZChkYXRhX2RpcltpXSwgIi8iLCAxMylbLDExXSkpKQogIHNhbXBsZV9uYW1lIDwtIHNvcnQoc2FtcGxlX25hbWUsIGRlY3JlYXNpbmcgPSBGQUxTRSkKICBwcmludChzYW1wbGVfbmFtZSkKICAKZm9yICh4IGluIHNlcV9hbG9uZyhzYW1wbGVfbmFtZSkgKSB7IAogIAogIHB5Y2xvbmVfaW5wdXRfc3Vic2V0IDwtIHB5Y2xvbmVfaW5wdXQgJT4lCiAgICBmaWx0ZXIoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCA9PSBzYW1wbGVfbmFtZVt4XSkgJT4lIAogICAgc2VsZWN0KC1jKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpKSAjJT4lIAogICAjIG11dGF0ZSh0dW1vcl9kZXNjcmlwdG9yID0gZmFjdG9yKHR1bW9yX2Rlc2NyaXB0b3IpLAogICAgIyAgICAgdHVtb3JfZGVzY3JpcHRvciA9IGZjdF9yZWxldmVsKHR1bW9yX2Rlc2NyaXB0b3IsIHRpbWVwb2ludHMpKSAlPiUgCiAgICAjYXJyYW5nZSh0dW1vcl9kZXNjcmlwdG9yLCBzYW1wbGVfaWQpCiAgCiAgIyBNYWtlIHRoaXMgcmVwcm9kdWNpYmxlCiAgc2V0LnNlZWQoMjAyMykKCiAgIyBEZWZpbmUgbGFiZWwgZm9yIHBsb3RzCiAgVGltZXBvaW50IDwtIGZhY3Rvcih4ID0gcHljbG9uZV9pbnB1dF9zdWJzZXQkdHVtb3JfZGVzY3JpcHRvciwgbGV2ZWxzID0gdGltZXBvaW50cykKICAKICAjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIwogICMgQ3JlYXRlIGJ4cCByZWZfY291bnRzCiAgICBwIDwtIHByaW50KGdncGxvdChweWNsb25lX2lucHV0X3N1YnNldCwgYWVzKHNhbXBsZV9pZCwgcmVmX2NvdW50cywgY29sb3IgPSBUaW1lcG9pbnQpKSArIAogICAgICAgICAgICAgICAgIGdlb21faml0dGVyKHdpZHRoID0gMC4xNSwgc2l6ZSA9IDAuNywgYWxwaGEgPSAwLjYpICsKICAgICAgICAgICAgICAgICBnZ3Bsb3QyOjpnZW9tX2JveHBsb3QoY29sb3IgPSAiYmxhY2siLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzaXplID0gMC4yNSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYWxwaGEgPSAwLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2VmID0gMCkgKyAjIHJlbW92ZSB3aGlza2VycwogICAgICAgICAgICAgICAgIHRoZW1lX1B1YmxpY2F0aW9uKCkgKyAKICAgICAgICAgICAgICAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gcGFsZXR0ZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYnJlYWtzID0gc29ydChuYW1lcyhwYWxldHRlKSkpICsKICAgICAgICAgICAgICAgICAjcm90YXRlKCkgKwogICAgICAgICAgICAgICAgIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTApKSArCiAgICAgICAgICAgICAgICAgc3RhdF9zdW1tYXJ5KGZ1bi55PW1lYW4sc2hhcGU9MSxjb2w9J2JsYWNrJyxnZW9tPSdwb2ludCcpICsKICAgICAgICAgICAgICAgICBsYWJzKHRpdGxlID0gc2FtcGxlX25hbWVbeF0sCiAgICAgICAgICAgICAgICAgICAgICB4ID0gInNhbXBsZV9pZCIsCiAgICAgICAgICAgICAgICAgICAgICB5ID0gInJlZl9jb3VudHMiLAogICAgICAgICAgICAgICAgICAgICAgY29sb3IgPSAiVGltZXBvaW50IikpCiAgICAKICAgICMgU2F2ZSB0aGUgcGxvdAogICAgZ2dzYXZlKGZpbGVuYW1lID0gcGFzdGUwKHNhbXBsZV9uYW1lW3hdLCAiLXJlZl9jb3VudHMucGRmIiksIAogICAgICAgICBwYXRoID0gcHljbG9uZV9wbG90c19kaXIsIAogICAgICAgICB3aWR0aCA9IDYsIAogICAgICAgICBoZWlnaHQgPSA1LCAKICAgICAgICAgZGV2aWNlID0gInBkZiIsIAogICAgICAgIHVzZURpbmdiYXRzID0gRkFMU0UpCiAgICAKCiAgIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMKICAjIENyZWF0ZSBieHAgYWx0X2NvdW50cwogICAgcCA8LSBwcmludChnZ3Bsb3QocHljbG9uZV9pbnB1dF9zdWJzZXQsIGFlcyhzYW1wbGVfaWQsIGFsdF9jb3VudHMsIGNvbG9yID0gVGltZXBvaW50KSkgKyAKICAgICAgICAgICAgICAgICBnZW9tX2ppdHRlcih3aWR0aCA9IDAuMTUsIHNpemUgPSAwLjcsIGFscGhhID0gMC42KSArCiAgICAgICAgICAgICAgICAgZ2dwbG90Mjo6Z2VvbV9ib3hwbG90KGNvbG9yID0gImJsYWNrIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2l6ZSA9IDAuMjUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGFscGhhID0gMCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29lZiA9IDApICsgIyByZW1vdmUgd2hpc2tlcnMKICAgICAgICAgICAgICAgICB0aGVtZV9QdWJsaWNhdGlvbigpICsgCiAgICAgICAgICAgICAgICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9IHBhbGV0dGUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGJyZWFrcyA9IHNvcnQobmFtZXMocGFsZXR0ZSkpKSArCiAgICAgICAgICAgICAgICAgI3JvdGF0ZSgpICsKICAgICAgICAgICAgICAgICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwKSkgKwogICAgICAgICAgICAgICAgIHN0YXRfc3VtbWFyeShmdW4ueT1tZWFuLHNoYXBlPTEsY29sPSdibGFjaycsZ2VvbT0ncG9pbnQnKSArCiAgICAgICAgICAgICAgICAgbGFicyh0aXRsZSA9IHNhbXBsZV9uYW1lW3hdLAogICAgICAgICAgICAgICAgICAgICAgeCA9ICJzYW1wbGVfaWQiLAogICAgICAgICAgICAgICAgICAgICAgeSA9ICJhbHRfY291bnRzIiwKICAgICAgICAgICAgICAgICAgICAgIGNvbG9yID0gIlRpbWVwb2ludCIpKQogICAgCiAgICAjIFNhdmUgdGhlIHBsb3QKICAgIGdnc2F2ZShmaWxlbmFtZSA9IHBhc3RlMChzYW1wbGVfbmFtZVt4XSwgIi1hbHRfY291bnRzLnBkZiIpLCAKICAgICAgICAgcGF0aCA9IHB5Y2xvbmVfcGxvdHNfZGlyLCAKICAgICAgICAgd2lkdGggPSA2LCAKICAgICAgICAgaGVpZ2h0ID0gNSwgCiAgICAgICAgIGRldmljZSA9ICJwZGYiLCAKICAgICAgICB1c2VEaW5nYmF0cyA9IEZBTFNFKQogIAogIAp9Cn0KCmBgYAoKCiMgVG90YWwgbnVtYmVyIG9mIG11dGF0aW9ucyBhY3Jvc3MgdGltZXBvaW50cyBhbmQgYmlvc3BlY2ltZW4gc2FtcGxlIHBlciBQYXRpZW50IGNhc2UKCmBgYHtyIGNyZWF0ZS1iYXJwbG90LXNhbXBsZSwgZmlnLndpZHRoID0gNiwgZmlnLmhlaWdodCA9IDUsIGZpZy5mdWxsd2lkdGggPSBUUlVFfQpmb3IgKGkgaW4gMTpsZW5ndGgoZGF0YV9kaXIpICkgeyAKICAKICAjIENyZWF0ZSBzYW1wbGVfbmFtZQogIHNhbXBsZV9uYW1lIDwtICB1bmlxdWUoYXMuY2hhcmFjdGVyKGdzdWIoIi5jYWxsLmNucyIsICIiLCBzdHJfc3BsaXRfZml4ZWQoZGF0YV9kaXJbaV0sICIvIiwgMTMpWywxMV0pKSkKICBzYW1wbGVfbmFtZSA8LSBzb3J0KHNhbXBsZV9uYW1lLCBkZWNyZWFzaW5nID0gRkFMU0UpCiAgcHJpbnQoc2FtcGxlX25hbWUpCiAgCmZvciAoeCBpbiBzZXFfYWxvbmcoc2FtcGxlX25hbWUpICkgeyAKICAKICBweWNsb25lX2lucHV0X3N1YnNldCA8LSBweWNsb25lX2lucHV0ICU+JQogICAgZmlsdGVyKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQgPT0gc2FtcGxlX25hbWVbeF0pICU+JSAKICAgIHNlbGVjdCgtYyhLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEKSkKICAKICAjIERlZmluZSBwYXJhbWV0ZXJzIGZvciBmdW5jdGlvbgogIHlsaW0gPSBtYXgocHljbG9uZV9pbnB1dF9zdWJzZXQkbXV0YXRpb25fY291bnQpCiAgCiAgCiAgICMgUmVuYW1lIGxlZ2VuZCBmb3IgdGltZXBvaW50cwogIFRpbWVwb2ludCA8LSBmYWN0b3IocHljbG9uZV9pbnB1dF9zdWJzZXQkdHVtb3JfZGVzY3JpcHRvcikKICAKICAjIFBsb3Qgc3RhY2tlZCBiYXJwbG90IAogIHByaW50KGdncGxvdChweWNsb25lX2lucHV0X3N1YnNldCwgYWVzKHggPSBzYW1wbGVfaWQsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeSA9IG11dGF0aW9uX2NvdW50LCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGZpbGwgPSBUaW1lcG9pbnQpKSArICAKICAgICAgICAgICAgICAgZ2VvbV9jb2wocG9zaXRpb24gPSBwb3NpdGlvbl9zdGFjayhyZXZlcnNlID0gVFJVRSkpICsKICAgICAgICAgICAgICAgZ2VvbV9iYXIoc3RhdCA9ICJpZGVudGl0eSIsIHdpZHRoID0gMC41KSArIAogICAgICAgICAgICAgICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBwYWxldHRlLCBicmVha3MgPSBzb3J0KG5hbWVzKHBhbGV0dGUpKSkgKyAKICAgICAgICAgICAgICAgdGhlbWVfUHVibGljYXRpb24oKSArIAogICAgICAgICAgICAgICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDg1LCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgaGp1c3QgPSAxLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgdmp1c3QgPSAxKSkgKyAKICAgICAgICAgICAgICAgbGFicyh0aXRsZSA9IHBhc3RlKHNhbXBsZV9uYW1lKSkgKyAKICAgICAgICAgICAgICAgbGFicyh4ID0gInNhbXBsZV9pZCIsIHkgPSAiVG90YWwgTXV0YXRpb25zIikgKwogICAgICAgICAgICAgICB5bGltKDAsIHlsaW0pKSAKICAKICAKfQp9CmBgYAoKYGBge3J9CnNlc3Npb25JbmZvKCkKYGBgCgoK